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ABSTRACT 

Zonally~averaged antocovariances of the five-day height tendency for 
Jags 0, 5, «+., 25 days are computed at 25 mb around latitudes 55N, 65K, 
and 75N during the period 5 October 195€ through 29 March 1959. Using 
the convolution theorem for each space-wave number, it is possible to 
determine the zonally=-averaged spectral density of the height-tendency 
field at each of these latitudes in terms of time waves having periods 
from 10 to 50 days. The results indicate that all tut about 10, of the 
variance of the height-changse field is accounted for by space waves one 
and two. For space wave number two, time waves number 3 and 4 (per 50 
days) give a significant contribution to the spectrum at the 905 confi- 
dence level. However, for space wave number one, time wave number one 
gives the most significant contribution to the spectrum, It is inferred 
from this, that cases of "explosive" warming, which mey occur as early 
as January ere correlated with intensification and/or movement associ- 
ated with space wave number one, at periods of approximately 100 days. 

The author wishes to express his sincere appreciation to Professor 
Frank L, Martin, Department of Meteorology, U. 5. Naval Postgraduate 


School, for his suggestion of the tonic and his encouragement in the 


preperation of this study. 
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phase angle of nth space wave at time t +U 
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raw spectral density value for time wave N 

smoothed spectral density value for time wave N 
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1. Introduction 

During the past ten years, meteorologists have become aware of the 
winter-time phenomenon known as "explosive" warming, which occurs in the 
Arctic and sub-Arctic stratosphere. During this same decade, the increased 
availability of higher-level soundings has made accessible to the meteor- 
ologist sufficient information to analyze the structure and causes of 
these unusual occurrences. 

The work of the McGill University group, primarily under Ff. K. Hare 
[1 960] and B. W. Boville [1960] has made it clear that there is a succes- 
Sion of cold troughs and warm ridges at the 25-mb level which eventually 
culminates in the development of an intense warm anticyclone in polar 
latitudes. Northern Canada is a preferred location for the occurrence of 
this phenomenon. This intense anticyclogenesis may be explained [Boville, 
1960| in terms of baroclinic instability in the polar night westerlies. 
Actually the development of the intense anticyclone represents a stage in 
the weakening of the polar night westerlies [Julian, 1959] which are re- 
placed in spring and summer by the easterly wind regime that finally en- 
velops the entire summer stretosphere. Boville and his associates [1961] 
have prepared an atlas of the 25-mb circvlation for the six-month pericd 
5 October 1958 through 29 March 1959. In these studies, Boville et al. 
list the Fourier analyses of the 25-mb contour field at five-day intervals 
for latitudes 55N, 65N, and 75N. For each of these latitude circles, they 


describe the Fourier analyses in the form 


_ « 
oe +2, cosinA-@,) (1) 


where the values of Cp, and ,, are listed in the atlas through space-wave 
number nine. Here % is the zonal mean contour height around the latitwle 
circle while Cy, and , are the amplituie and phase angle, respectively, 
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of the nth harmonic (n=1,.....,9). 

Spectral analysis has become a useful diagnostic tool of meteorolo- 
gists, primarily for the purpose of detecting preferred periodicities in 
the eddies or perturbations under consideration. The number and variety 
of such stucies in certain areas of statistical meteorology are too great 
to attempt to describe them here. Suffice it to say that the basic methods 
were largely developed by J. i. Tukey, and a concise swmary of their ap- 
plications to Meteorology can be found in the excellent textbook by 
Panofsky and Brier fiose]. 

Actually, Boville [1960] has vcublished the spectral analysis of the 
25-mb contour height field at latitude 60N for the month of October 1958. 
In that analysis, one of the conclusions reached was that most of the 
variance in the contour height field at 25 mb is explained by space-wave 
numbers one and two. A similar conclusion is reached in this paper in 
regard to the variance of the height tendency. 

The objective of the present paper may be understood from the local- 


time differentiation of equation (1) 


Z “2 IC cos(nA-d,| + Cadnsin(N—d,) (2) 


n=4 

where 32¥ =O z/ot _ 07/0 5 . From equation (2), it may be 
seen that time-changes in the perturbation field may be ascribed to terms 
involving cs the local-time change in the harmonic wave amplituce, and 
dp » the local-time change of the phase angle. These two types of terms 
in equation (2) may be described as the intensification and movement con- 
tributions, respectively, for the nth wave. A time-spectral analysis of 
dz/dt, zonally averaged with respect to longitude, was then performed 
using the five-day rates of change Ce and d, . (For the mathematical de- 
velopment see section 2.) An advantage of the spectral analysis of the 
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height-tendency field over that of the height field alone is that the 


former can indicate periodicities due to intensification and to movement 





of the harmonic waves, whereas these connections are not so clear in the 


latter case. 





2. Mathematical development 
The autocovariance of z* atlases, 1S gived [Blacknan and Tukey, 


1958] by the formula 
, We 
R(t) =Gam 2 Aare 26 st) dt (3) 


-T/2 
where #™* is the zonal height perturbation, and includes the contribution 


from all space waves around the latitude circle, T is the averaging per- 
iod over which the autocovariance is determined, T is the lag, and t is 
an interval of time, also measured in five-day units. Using the truncated 
time series available for this study, the time integration may be replaced 
by a summation over NK equally-spaced time intervals, at lag T= KAt. 


Then the integral (3) may be rewritten in the form 
ar N-k 


Mik) = 1 i) [2% ZR(taT) | SA. Omi Oa 
* t=1 


where R(K) is the zonally-averaged time autocovariance. In equation a). 
A represents the longitude at any point within a latitude circle, Here, K 
is tne number cf lags in terms of intervals of At =35 days, J] is the num- 
ber of five-day height changes available in the time series, and M is the 
maximum number of time lags used in this study. Obtaining the five-day 
perturbation height tendency z2~ ( t+U) at time [+T dy a equation 


analogous to (2), we may write, 
9 


Ze") =) (é, cos (NA- W,) + G. U, sin (NA -),)} (5) 
nel 


where @, and W, are analogous to Cn and o, of equation (1), Hence, 
from equations (2), (4), and (5), the zonally-averaged time autocovariance 


for space wave n has the form 





. — Pa P fy 
Friis i [ont Ge @ tf on (Y,-O,) + 
f xm | | 
| ; 


[aCe = Call] sin(U,-O)} (6) 


In equation (6) the subscript n indicates that R, (kK) is independent of all 
space waves other than the nth. In the derivation of (6), considerable 
simplification has been made possible as a result of orthogonality prop- 
erties, that is by deleting terms which vanish upon integration around a 
complete latitude circle. A proof of equation (6) is given Appendix I. 
Since discrete data and a finite time series have been used, the maximum 
lag Mat linits the maximum resolvable period to be [,,.,=2 MAt 
which may be called the fundamental period. In this study the maximum 
number of lags has been taken as M = 5 so that the fundamental period will 
be 50 days. If K is allowed to range over the values 0, 1, 2, «woe, M, 
various values of Ry {K) can be generated for the different space-wave 
numbers, Blackman and Tukey [1958] recommend, when using a truncated 
time series, that the maximum number of lags M should not exceed approx- 
imately 10% of the length of the record. Since the total available number 
of five-day mean values of dz We ot was 35, the maximum number of lags 
was taken to be five. This lag number was considered to be in reasonable 
accord with the recommendation of Blackman and Tukey, and yet sufficiently 
large to give resolution for the long-period waves which occur in this 
meteorological regime. | 

The definition of the spectral density at time-wave number N (Kahn, 
1957] is the average contribution to the variance from a unit wave number 
interval centered at wave number N. Since the maximum period under con- 


sideration is 50 days, the concept of a time-wave number N is to be inter- 





preted in terms of frequencies as follows: 

1 wave per 50 days for N= 1, 

N waves per 50 days for wave number N, N = 5, 
Thus the frequency - may be expressed as 

{y= N : N<é Mi (7) 
Z2MAt 

and spans the frequency-interval (N -3+)/2MAt to (N+3)/2MAT. The 
corresponding period for this time-wave is the reciprocal of fue The 
highest frequency fe resolvable when using such discrete data at finite 


intervals occurs for N = M and is 
fe 3a 
F ZS 
f. is often called the folding or Nyquist frequency, and in the present 
study is one cycle (10 days i ; 


Panofsky and Brier [1958] give the following equation for raw spec- 


tral density estimates for time-wave number N: 


1 
Si = 2{4[Ro+ Cr Rev) + [Rew COS TNA}, 10 =™M-1. (8) 
Mines L M 
In practice, 5, is generally considered to be distributed uniformly over 
the unit wave-number interval N - 3 to N+ <4, except at N= 0 and N=™M, 
For the special cases N = O and N=M, the corresponding estimates — 
S, and S,, follow from equation (8) multipled by one-half because only 
one-half of a wave=number interval is available at N= 0 and N=M, 
These unsmoothed raw spectral estimates may be smoothed by "hanning" 


{ Blackman and Tukey, 1958], in the following manner 


Dy =e NH ee Oa (9) 
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except at the extreme wave numbers N = O and N = M where the average is 


modified to 


55-5. 5,405 Sts ees, +-5 Sm, (10) 


These smoothed estimates may be considered to be centered at wave number 
N, or frequency $y given by equation (7). 

Note that the integrated value of all the spectral estimates for a 
particular space-wave n gives the total time variance, for this wave. I 


other words, by definition, the relationship 


= 


vo) 


2 
Oh, =) SypAN 4) 


N= 
holds, although carrying the integration only to N = 5 neglects some of 


the higher frequency spectral contributions. In addition, "aliasing" or 
folding of high frequency energy back upon the Nyquist frequency has been 
neglected. However the spectral densities at higher space-wave numbers 
(table 2) suggest that this is not important in the present problem. 

This has been inferred since the contributions to the time variance for 
these wave numbers appears to act as "noise", with no preference for the 


higher time-wave numbers. 





3. Computational procedures and interpretations 

As mentioned in section 2, the maximum number of time lags chosen 
was M= 5. This choice gives a resolution of the spectrum of dz/odt 
to a maximum frequency of 1/10 cycle (aay). 

All computational processes were performed on the CDC~1604 computer. 
All programs devised for this study were written by the author or modified 
from those available in the U. S, Naval Postgraduate School computer cen= 
ter library. The three latitudes chosen for analysis were 55N, 65N, and 
75N. Table 1, a, b, c gives the zonally-averaged autocovariances for each 
space-wave number computed by equation (6) at each of the three latitude 
circles, The smoothed density spectra for each space-wave number were 
then computed according to equations (8), (9), and (10) and the results 
displayed in table 2, a, b, ¢. 

Note that in table 2, a, b, c periodicities are presented down the 
left-hand column, while the space-wave numbers n= 1, ..., 9 appear along 
the top of the table. From the rows denoted "P, R. Var.", it may be seen 
that the space waves n = 1 and n = 2 together account for approximately 
90% of the variance of the zonally~averaged five-~lay height changes. 
While the percent reduction of the variance of 0z/dt accounted for by 
space wave 2 exceeds that of wave 1, the contribution of the latter in- 
creases as one proceeds from latitude 55N to 75N. The contribution by 
waves 3 and 4 accounts for less than 10% at latitwies 55N and 65N, and 
for less than 5% at latitude 75N. Thus, in this study, only the time- 
spectra of space waves one and two have been analyzed for significant 
periodicities. 

For space waves n= 1 and n= 2, fiducial confidence limits at the 
10% and 90% significance level have been taken from Blackman and Tukey 


[1958]. These limits are based upon the VA distribution, where 
8 





LAG 
DAYS 
4 
0 4218 
9 8.86 
10 243 
15 =3 37 
20 «86 
a, -3..46 
0 80.78 
2 23 62 
10 “1.55 
ES ~4,.95 
20 1.33 
a (S95) 
O® 103067 
2 28 .08 
10 7.89 
1. el .79 
20, 4.05 
25 11.88 
Table 1. 
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Oe 
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Zonally averaged autocovariances of 02/dt 
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(a) 55N 


FREQUENCY SPACE WAVE NUMBERS 
CYCinS piu 
1 2 e 4 5 6 a g 9 
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42.0 5300 


2.6 


1.1 0.5 


Zonally-averaged spectral densities 


0.1 0.1 0.1 


for space-wave number n=1, 


20059, and wave frequency N=0,1/50,...,5/50 cvcles/day. 


sec- 


tions (a), (b), and (c) are for latitudes 55N, 65N, and 75N. 


The units are in terms of 10*(feet/5 days)*. 
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fe @ = MA\/™ is the number of degrees of freedom. In this in- 


vestigation we have oa = 35 and M = 5, so that f A714. For these values, 
the 10% and 90% fiducial limits interpolated from those given by Blackman 
and Tukey nese. p.208] are 0,55 and 1.51 respectively. 

Ficure 1 illustrates the time-spectral analysis at 65N, for space 
waves n = 1 and n = 2 (the latter show by the solid line). From this 
diagram it may be seen that for n = 2, the mean variance averaged over all 
periods is 60.1 x 10*(e4/5 days )*. The fiducial limits have been con- 
structed relative to this mean. Thus it is evident that for n = 2, waves 
of period 10 days and 12.5 days give significantly large contributions at 
a 90% confidence level. On the other hand, at 65N space wave 2 at a per- 
iod of 50 days, gives significantly small contributions to the variance. 

The same type of analysis for n = 1 at latitude 65N shows that the 
mean variance for space wave n = 1, averaged over all periods is 13.5, in 
units of 104(£4/5 days )*, so that a time-wave of period 50 days gives a 
significant contribution at a 90% confidence level. It should be recalled 
that a period of 50 days in figure 1 actually spans a range of periods 
from 33.3 to 100 days, and such waves could include those associated with 
explosive warming. This is so, since strictly speaking, the largest value 
of dC, /ot which contributes to the variance of 0Z/ at could occur at 
a time which is either three-fourths or one and three-fourths of a period 
after the onset of the autumnal cooling season in polar latitudes. 

In table 3 the mean spectral densities at latitudes 55N and 75N, 
averaged over all frequencies in table 2, are given for space waves n = 1 


eon = 2e 
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2 
Table 3. Mean spectral densities in units of 10*(feet/5 days) at lati- 
tudes 55N and 75N for space-waves n= 1 and n= 2. 


Using the fiducial confidence limits for the es distribution pre- 
viously mentioned (0.55 and 1.51, respectively), it is apparent from table 
2 a; c} that there exists no significant preferred peak in the spectral 
densities at the 90% limit for either latitude 55N or 75N, 

For space wave n= 1, at each of the latitude circles of table 2, 
there is a steady decrease of spectral density, S,, with increasing time- 
wave N=1. 

For the case of isotropic turbulence, MacCready [1953] has tested the 


applicability of a spectral function of form 
= 
oe CaN 


where N is the time-wave number. This test is essentially based upon a 
theory of isotropic turbulence set forth by Heisenberg [1948] in which 
p= 5/3. Forn=1, and using the results of table 2, computation of the 
values of the exponent p gave the following: 

po, =p = 0857 

Oe p= 0.91 

Woe, «op = O.64 
Chiu [1960] has obtained somewhat similar results for the exponent p using 
the spectrum of the zonal and meridional components of the wind at Belmar 
and Cocoa, and finds values of p of the order of 1. 

On the other hand, the results of this study for space wave n = 2 do 
1s 





not show the decrease of spectral density with increasing time-wave num—- 
ber. Hence for n = 2, the peak spectral density which occurs in each case 
of table 2 at either time-wave 3 or 4 is even more likely to be signifi- 
cant of the average state of synoptic events. 

Space-wave n = 1, with its maximum spectral contribution at time-wave 
N= 1 is capable of explaining those perturbations of longer period asso- 


ciated with "explosive" warming. 
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4. Conclusions 

(1) Space waves 1 and 2 contribute 90% of the variance of the height- 
tendency field at latitudes 55N, 65N, and 75N, 

(2) The chief contributors to the spectrum of szfat occur for 
time~vavesN = 3 or N= 4 (periods 12.5 and 16.7 days), which have the same 
general periodicity as the average perturbation period at a fixed station 
(see, for example, Austin and Krawitz, 1956]. 

(3) The chief contributor to the time variance for space wave n = 1 
occurs for time wave N = 1, or vossibly N = -~5, which corresponds to periods 
of the order of 50 days or longer. Such veriods may be sufficiently long 
to account for the "explosive" warming which sometimes occurs as early as 
January. 

(4) For space weve n = 1, the time spectral density was distributed 


with time wave number according to 


=O 
Oyo N 


which compares reasonably with the results obtained in other spectral 


studies [c.f. Chiu, 1960]. 
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APYENDIX I 
DERIVATION OF THE ZONALLY-AVERAGED TIME AUTOCOVARIANCE 
The equation for the zonally-cveraeged time autocovariance of the 
height-tendency field at lag ‘i may be expressed in the form shown in 


equation (4) 


an ‘R-K 


A (kK) = zt) Zz (tat)]SA, ww 


where 7] is the maximum number of discrete values in the sample, and Ras) 
is the zonally-averaged time autocovariance. 
When equation (4) is expanded using equations (2) and (5), and upon 


interchanging the order of summation, equation (12) is obtained 
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where the orthogonality property has also been used. Equation (12) repre- 
sents a double integral; first, with respect to longituwle, second, with 
respect to time. The integral with respect to longitude gives the zonal 
average. On integrating around the latitude circle from 0 to 2TT, all 
cosine and sine terms involving the argument (2 nA — Dn — Wn} 


integrate to zero. Therefore, the result obtained is 
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Note that R n(K) is the zonally-averaged autocovariance at a time-lag 
of KAT days, for the nth space-wave; no other space-wave appears in 
the expression for R n(K) » the interval between successive date is 


such thet Al= 5 days. 
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